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Excitation Energies in Density-Functional Theory 
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Specific matrix elements of excliange and correlation kernels in time-dependent density-functional 
theory are computed. The knowledge of these matrix elements not only constraints approximate 
time-dependent functionals, but also allows to link different practical approaches to excited states, 
either based on density-functional theory, or on many-body perturbation theory, despite the approx- 
imations that have been performed to derive them. 
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Density- functional theory (DFT), as proposed by Ho- 
henberg, Kohn, and Sham has been highly success- 
ful in the analysis of interacting-particle ground states. 
Many efforts have also been devoted to DFT- 

based schemes for excited states. The computation of ex- 
citation energies from DFT total-energy differences 
a procedure known as ASCF, is formally restricted to 
transitions between the lowest states of different symme- 
tries. The ensemble-density- functional theory ||^, a more 
general path towards excitation energies, has not been 
used for actual calculations, to our knowledge. Recently, 
the feasability of excitation-energy computations relying 
on two other, widely applicable, DFT-based formalisms 
has been demonstrated for atoms and small molecules. 
The first Q ^ starts from the extension of DFT to time- 
dependent (TD) phenomena (TDDFT). The sec- 
ond [p^|-|l2| , due to Gorling and Levy, builds a perturba- 
tion theory (GLPT) in the difference between the many- 
body and the second-quantized Kohn-Sham (KS) Hamil- 
tonians, where the parameter of the perturbation is the 
coupling constant of the particle interaction, in such a 
way that at each order, the exact density is recovered. 

As an alternative to these DFT-based efforts, one may 
start from many-body perturbation theory and perform 
partial resummations of diagrams such as to build a 
screened interaction between dressed particles (quasi- 
particles) At the lowest order in the screened in- 

teraction, one obtains Hedin's GW approximation jl^] 
to quasi-particle energies (one particle is added or sub- 
tracted to the system), while the energy of excited states 
for which the number of particles is conserved can be 
deduced, in a subsequent step, from a Bethe-Salpeter 
(BS) equation describing the interaction between quasi- 
particles (e.g. electron- hole pairs). The application of 
these techniques to real materials is very demanding . 

In the present Letter we explore the relationships be- 
tween excitation energies derived from time-dependent 
density-functional theory, the Gorling-Levy perturbation 
theory, and the screened-interaction many-body pertur- 
bation theory. If excitation energies were derived with- 
out approximations, in the three formalisms, the results 
should be identical. Because of the specific theoretical 
developments, the most obvious simplifications are dif- 



ferent, so that the practical schemes derived from these 
formalisms also differ. We find that some approximations 
used for practical calculations leave a connection between 
the approaches, at variance with the usual adiabatic 
local-density approximation (ALDA) in TDDFT [|-§, 
that leads to a different physical picture. 

The expressions from different approaches are linked 
thanks to a new technique for computing selected ele- 
ments of TD functional kernels at resonance, that is, at 
the frequencies corresponding to differences in KS eigen- 
values, for which the independent-particle susceptibility 
of the KS system is resonant. It is first applied to the ex- 
act exchange kernel, whose matrix elements appear in a 
simplified TDDFT treatment of excitation energies based 
on a Laurent expansion, and are found identical to the 
first-order corrections to KS eigenvalues differences, in 
the GLPT. The knowledge of these matrix elements im- 
poses a new constraint on approximate functionals. The 
same technique is then applied to an explicit XC func- 
tional that includes an approximate correlation contri- 
bution. It is shown that the corresponding Laurent- 
expansion TDDFT excitation energies can be split ac- 
cording to powers of the screened interaction, as in the 
GW^-BS approach, and that they contain the GW quasi- 
particle eigenvalue shifts. 

For sake of simplicity, we will consider systems made 
of spinless particles (or equivalently, fully spin-polarized 
systems) and focus on finite systems for which the energy 
levels of interest are non-degenerate (KS wavefunctions 
are real) . The introduction of spin dependence and com- 
plex wavefunctions is technically straightforward, and 
does not modify the conclusions reached here. We use 
Hartree atomic units throughout. For integrals of the 
form / gi(ri)/(ri,r2)52(r2)dridr2, in which gi and 32 
are not wavefunctions, but densities or products of wave- 
functions, we will use the notation {gi\f\g2}- 

In TDDFT, the lowest-order Laurent expansion of the 
XC kernel and KS susceptibility, in the difference be- 
tween the actual excitation frequency, fi, and the en- 
ergy difference between corresponding KS eigenenergies, 
uj^^ = e^^ — e^^, allows to deduce approximate excita- 
tion energies as ^ : 
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The double index v = (j, k) labels a single-particle 
transition from the occupied level k to the unoccu- 
pied level j, $i/(r) is the product of KS wavefunc- 
tions (/>fc(r)(/)j(r), Vc is the Coulomb interaction kernel 
( | ri-r2| ) ^^'^ /xc(ri,r2;Ci;) is the frequency-dependent 
XC kernel. 

In Eq.(|l|), the diagonal matrix elements of the 
Coulomb and XC kernel for the transition between 
occupied and unoccupied states correct the difference 
in KS eigenvalues. The Coulomb kernel gives an 
exchange-like interaction between the KS electron and 
hole [</'fc(ri)0,(ri)^;:^</>fe(r2)0j(r2)]. The XC kernel 
is easy to evaluate in the ALDA (that ignore the func- 
tional dependence of the XC potential Vxc on densities at 
different times as well as different locations): It induces 
only a local, attractive, static electron-hole interaction 
[0fc(ri),/),(ri)5(ri ~ r2)^(ri)0fc(r2)</),(r2)]. Then, as 
detailed by Cassida [Q, quasi-particle energies (i.e. exci- 
tation energies, to which quasi-particle interactions make 
no contribution, like ionisation energies) are identified 
with unmodified KS energies [|9|. 

Unlike the ALDA, we evaluate exactly the exchange 
contribution. Since we will use later the same technique 
for the correlation contribution, we treat them now on 
the same footing. The Keldysh formalism, as applied re- 
cently to the TDDFT case by van Leeuwen (we follow 
his notations), leads to the following explicit expression 
for the TD kernel, second derivative of the action A^c, in 
the temporal (pseudo-time) domain: 
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In this equation u and Uxc, respectively, denotes the time- 
domain KS potential and XC potential. 

This general expression can be specialized to the TD 
linear response of an otherwise static system, in the fre- 
quency domain, and inserted in the integral needed in 
Eq.(0): 



{$,|/xc(c^)|$.} = {$. 
/ S^Ax,{lj = 0) 



dr5Uxc(r5) 
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Note the appearance of the (static) XC potential Vxc (r) . 
(5ws(r;w) is a change of KS potential in the frequency 
domain, xo(ri,r2;'^) is the independent-particle suscep- 
tibility of the KS system, and Xo^ is the effective inverse 
of xoi defined on a restricted space excluding constant 
functions M . The spatial dependence of potentials and 



susceptibilities is implicit in Eq.(||), as well as most inte- 
grals over spatial arguments. A^duJ = 0) and n(r; lo — 0) 
are the time-independent part of the action integrand and 
the density, respectively M. 

When evaluated at lo — lJ^^, Eq.(|3|) gives the XC ker- 
nel contribution to in Eq.(|l|). This requires the in- 
version of xo around the frequency w^^, and its sub- 
sequent application to ^,^(r). The explicit expression 
of Xo is well-known (see Eq.(14) of Ref. |^), and ex- 
hibits a resonance at that frequency. The correspond- 
ing pole, can be singled out before performing the inver- 
sion: xo(ri,r2;t^) = '^"t-^^s"''^ + Xi.(ri,i'2;^), where 
the non-resonant contribution Xu is implicitly defined by 
this equation. The Sherman- Morrison formula pl[ allows 
the analytical inversion of such a decomposition. The ap- 
plication of Xo ^ to <i>y, denoted v^, gives: 

v,{uj) = Xo\^)\'^.} 
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A similar expression involving also uj — uj^^ can be de- 
rived in the case of v^{—uj). The presence of the fac- 
tor {uj — uj^^) emphasizes the anti-resonant behaviour 
of Xo^- This characteristics of Eq.(|^) allows to extract 
the leading order of the frequency expansion of Eq.(|^) 
around lu^^, for XC actions whose explicit expression 
is known as a function of the KS potential through KS 
wavefunctions. Indeed, after technically straightforward 
calculations within perturbation theory (up to second- 
order wavefunctions), one finds 



(S^n(r; lu — 0) 
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and, using the explicit exchange action of Ref. fl^. 
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where is the Fock non-local operator evaluated with 
KS wavefunctions and nk{v) is the square of KS wave- 
function 4)k{'r). Combining Eqs.(j|),(||) and (^), we de- 
rive the following matrix element of the exchange kernel 
at resonance (in the limit of w — > i^^^), valid for all pairs 
(j, k) of unoccupied and occupied levels associated with 
a composite index v: 



{$.i/x(c.:^^)i$.} = 
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- {nk\Vc\nj} . 

This explicit expression for selected matrix elements of 
the TD exchange kernel provides, as such, a new con- 
straint on approximate TD exchange functionals. Note 
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the drastic simplification from the general matrix ele- 
ments of the TD exchange kernel derived by Gorling . 
The combination of Eqs.(|^) and (Q), neglecting the cor- 
relation contribution, gives 



n 
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This expression is identical to the one obtained in the 
first-order GLPT (see Eq.(A3) of Ref. and also 

Ref. p!c|]): Excitation energies from the TDDFT with ex- 
act treatment of the exchange kernel, in the Laurent ap- 
proximation, and from first- order GLPT are equal |^^ . 

The KS eigenenergies e*^^ and e^^ are expectation val- 
ues of the KS Hamiltonian, that includes the exchange 
potential. Hence the two expectation values of the differ- 
ence between the Fock operator and the exchange poten- 
tial in Eq.(||) induce a shift of KS energies and e^^, 
present in u)^^ , in the direction of the Hartree-Fock quasi- 
particle energies. So, unlike in the ALDA approach, the 
quasi-particle energies from TDDFT with exact treatment 
of the exchange kernel are shifted from the KS eigenener- 
gies. The fourth term in Eq.(^) describes an unscreened 
Coulomb electron-hole attraction, to be contrasted with 
the effective interaction present in ALDA. However, the 
third term of Eq. (H) is the same exchange- like interaction 
as in the ALDA case, present in Eq.(|l]). 

The dynamical treatment of correlation is more diffi- 
cult than that of exchange. We will not pursue the inves- 
tigation of GLPT, that would involve higher-order terms. 
However, we will establish a link between TDDFT and 
the screened-interaction many-body perturbation theory. 
At variance with the decoupling between GW and BS 
procedures, TDDFT computes the excitation energies in 
one operation. This might be a crucial advantage for 
practical applications. Our aim will be to show, at the 
adequate approximation levels, that TDDFT contains 
the ingredients characteristics of the GW approach, and 
to point out how it accounts for the additional BS terms. 

The GW quasi-particle shift, evaluated in a non- 
self-consistent, diagonal approximation, based on KS 
eigenenergy and eigenfunctions, leads to the following ex- 
citation energies: 



.fe|S«^(eP)-t;.e|0fc) 



(9) 



Y^GW jg ^YiQ self-energy operator in the GW approxima- 
tion: 

sG^(ri,r2;e)=y' ^e-^'UG{r,,r2; e - e')W{rur2; e'), 

where the Green's function G is built from KS wavefunc- 
tions, and W is the dynamically screened Coulomb inter- 
action {W = £^pjiVc, with £rpa = 1 — VcXo being the 



RPA dielectric screening). The e'-integral can be further 
transformed to a contour integral in the complex plane. 

In Eq.(^, the difference between KS eigenenergies is 
shifted by expectation values of the difference between 
the self-energy operator and the XC potential. The ne- 
glected BS corrections would make quasiparticles inter- 
act through an unscreened exchange-like term (compare 
with Eq.(|l|)) as well as through a screened Coulomb in- 
teraction. Excitonic effects would be described at that 
level. The screened Coulomb interaction is to be con- 
trasted with the effective, local, electron-hole interaction 
of ALDA-TDDFT, and the unscreened Coulomb interac- 
tion of first-order GLPT. 

A formally exact correlation functional can be ob- 
tained by the combination of the adiabatic connection 
method (based on a coupling constant integral) and the 
fluctuation-dissipation theorem . In the context of the 
DFT computation of van der Waals energy, Dobson [2^ ] 
proposed a simplified expression, 

^"^Jo'^^Jo ^^'[^c(Xo(«")-XRPAa(«u))] ,(10) 

in which the true susceptibility present in the exact 
expression is replaced by the RPA one, although the in- 
tegral over the coupling constant A is retained. 

The RPA susceptibility is a simple functional of xa 
[ Xrpa,a = Xo(l - AVcXo)"^ ]• This latter quantity is a 
functional of the KS potential, itself a functional of the 
density. We can thus compute the corresponding corre- 
lation potential and kernel. Thanks to an integration by 
parts, we obtain that the correlation potential does not 
exhibit the integration over the coupling constant A: 
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where the same screened interaction W as in the GW 
approximation appear. For the computation of excited 
states, we generalize the correlation functional Eq.(p^) to 
a TD action, and, neglecting a term in the second power 
of the screened interaction, derive 
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For comparison with the GW approximation, it is crucial 
that the coupling constant is absent, as in Eq.(pT|). The 
second derivative of xo with respect to TD changes of 
potential can also be evaluated (perturbation theory up 
to third-order wavef unctions), close to the resonance: 
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where the same Green's function G as in the GW ap- 
proximation appears. The combination of Eqs.(12) and 
(H), needed in Eq.(|), gives the difference between an 
unscreened GVc term and a GW term. We introduce 
this result in Eq.(||), and find that the main contribu- 
tion of the first term (GVc ) cancels the unscreened Fock 
operator that shifts the KS eigenvalues in Eq.(||). 

We can now discuss the link between GW-BS excita- 
tion energies and TDDFT excitation energies obtained 
from an exact exchange kernel and an approximate cor- 
relation kernel : i) as expected, the crucial GW terms 
are contained in Eq.(^; ii) the inclusion of correlation 
in TDDFT transforms the unscreened single-orbital ma- 
trix elements of Eq.(|) to dynamically screened matrix 
elements, as in Eq.(|£); iii) the subtraction of the same 
matrix elements of v^c appear in both schemes, compare 
Eq.(|[^) to Eq.M; iv) however the contour integral de- 



duced from Eq.( 12) is shifted from the contour integral in 
the GW approximation, and includes contributions from 
additional poles ; v) a term in the second power of the 
screened interaction has been neglected from Eq.( |lO|) t o 
Eq.(|l2|); vi) the RPA susceptibility was used in Eq.(|lO[). 
Points iv)-vi) connects to BS as well as 0{W") (with 
n > 2) corrections to Eq.(||). Understanding point iv) 
also requires a careful analysis of our TD correlation ac- 
tion (especially the contour integral to be used with the 
fluctuation-dissipation theorem). A more detailed discus- 
sion of these aspects will be left for a future pubhcation. 

In this paper, we have examined the time-dependent 
kernels deduced from the exact exchange functional and 
an approximate correlation functional, at resonance. 
Links between approximate expressions for excitation en- 
ergies from time-dependent DFT, the Gorling-Levy per- 
turbation theory, and Hedin's GW approximation have 
been found: The difference between Kohn-Sham eigen- 
values are corrected by one-particle shifts, followed by 
(unscreened or screened) interaction effects. We also 
noted that in the adiabatic local density approximation 
to time-dependent DFT, no one-particle shift is present, 
and the link with the two other techniques is lost. 
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